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Magnetic and electronic properties of the Hubbard model on the Bethe and fee lattices in infinite dimen- 
sions have been investigated numerically on the basis of the dynamical coherent potential approximation 
(CPA) theory combined with the harmonic approximation (HA) in order to clarify the quantitative as- 
pects of the theory. It is shown that the dynamical CPA+HA reproduces well the sublattice magnetization, 
the magnetizations, susceptibilities, and the Neel temperatures (Tn) as well as the Curie temperatures 
calculated by the Quantum Monte-Carlo (QMC) method. The critical Coulomb interactions (Uc) for the 
metal-insulator (MI) transition are also shown to agree with the QMC results above Tn. Below Tn, Uc 
deviate from the QMC values by about 30% at low temperature regime. These results indicate that the 
dynamical CPA+HA is applicable to the quantitative description of the magnetic properties in high dimen- 
sional systems, but one needs to take into account higher-order dynamical corrections in order to describe 
the MI transition quantitatively at low temperatures. 

KEYWORDS: Dynamical CPA, Curie temperature, Nee! temperature. Effective Bohr magneton number. Metal- 
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1. Introduction 

The magnetism of correlated electron systems with in- 
termediate Coulomb interaction strength has been one of 
the long standing problems in the solid-state physics be- 
cause these systems show the complex properties known 
as the itinerant vs. localized behaviors.^' The ground- 
state magnetizations of Fe, Co, Ni, for example, show the 
noninteger values in unit of the Bohr magneton number 
being characteristic of a band model, while their finite- 
temperature properties such as the magnetization vs. 
temperature curves and the Curie- Weiss susceptibilities 
are well explained by a Heisenberg-type localized model. 
The Hartree-Fock theory and simple perturbation the- 
ories could not explain their magnetic properties since 
these theories overestimate the magnetic energy at finite 
temperatures. 

Because of the difficulty in describing the itinerant vs. 
localized behavior in the magnetism of transition met- 
als, interpolation theories have been developed by many 
investigators. Cyrot'^^ adopted the functional integral 
method'*-' to the Hubbard model and derived the P-T 
phase diagram of metal-insulator (MI) phase transition 
at finite temperatures. Hubbard^-* and Hasegawa^-' inde- 
pendently developed the single-site spin fluctuation the- 
ory (SSF) on the basis of the functional integral method. 
The theory transforms the electron-electron interaction 
into a time-dependent random potential by introducing 
the time-dependent charge and exchange random ficti- 
tious field variables, and treats the potentials using the 
coherent potential approximation (CPA).'') The theory 
interpolates between the weak and strong Coulomb in- 
teraction limits, and explained qualitatively the localized 
vs. itinerant behaviors of magnetism in transition metals. 

The SSF is based on the static approximation (SA) 
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which neglects the time dependence of the fictitious 
fields. Though the SA is exact in the high tempera- 
ture limit, it reduces to the Hartree-Fock theory at the 
ground state. Thus it does not take into account the elec- 
tron correlations as found at the zero temperature by 
Gutzwiller,^' Hubbard,^) and Kanamori.-'^^^ Kakehashi 
and Fulde*"'^) proposed a variational theory at finite tem- 
peratures which takes into account the ground-state elec- 
tron correlations, and showed that the correlations at 
finite temperatures can reduce the Curie temperatures 
obtained by the SSF by a factor of two. 

Towards more quantitative theory, Kakehashi^^^ pro- 
posed the dynamical CPA which completely takes into 
account electron correlations at finite temperatures 
within the single-site approximation, and clarified the 
basic properties of the theory using a Monte-Carlo tech- 
nique. Later he proposed an analytic theory of the dy- 
namical CPA combined with the harmonic approxima- 
tion (HA),^'^'' and examined the dynamical effects on the 
itinerant ferromagnetism. In the HA,*^) we neglect the 
mode-mode couplings between the dynamical potentials, 
and take into account the contributions from the dynam- 
ical potentials with the same frequency independently. 
The HA is known to be exact up to the second order 
in the weak Coulomb interaction limit, and is known 
to describe quantitatively the Kondo limit in the strong 
Coulomb interaction regime. The dynamical CPA has 
recently been proved^' to be equivalent to the many- 
body CPA in the disordered problem,*^) the dynamical 
mean-field theory in the metal-insulator transition in in- 
finite dimensions, and the projection-operator method 
CPA for excitations in solids. 

More recently, we proposed the first-principles dy- 
namical CPA^°) on the basis of the tight-binding linear 
muflintin orbital (TB-LMTO) Hamihonian2i'22) with 
the LDA (Local Density Approximation) +L/ type in- 
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traatomic Coulomb interactions, ^'^^ and clarified the fer- 
romagnetism of Fe, Co, and Ni,^^'^^' as well as the sys- 
tematic change of excitation spectra in 3d transition 
metal series. ^^■^'^^ 

Although the numerical results of the first-principles 
dynamical CPA combined with the HA seem to be rea- 
sonable and explained well a systematic change of the 
XPS and BIS data,^^) the validity of the dynamical 
CPA-I-HA has not yet been examined in details. The pur- 
pose of the present paper is to clarify the quantitative 
aspects of the dynamical CPA+HA from the numerical 
point of view. For this purpose, we have calculated the 
antiferromagnctic properties of the half-filled Hubbard 
model on the Bethc lattice in infinite dimensions and 
the ferromagnetic properties of the non half-filled Hub- 
bard model on the fee lattice in infinite dimensions'**^ 
on the basis of the dynamical CPA+ HA. We clarify the 
validity of the theory by comparing the results of calcu- 
lations with those of the Quantum Monte-Carlo (QMC) 
method.'®"^^^ We will conclude that the dynamical CPA 
approach is useful for the quantitative understanding 
of the magnetic properties of solids in high dimensions 
where the single-site approximation works well. 

Outline of the present paper is as follows. In §2, we 
summarize the dynamical CPA combined with the HA 
and extend the theory to the antiferromagnctic case with 
the bipartite lattice. In §3.1 we present the results of 
calculations for the antiferromagnctic states in the half- 
filled Hubbard model on the Bethe lattice. It is shown 
that the zeroth approximation to the dynamical CPA, 

1. e., the SA quantitatively describes the Neel temper- 
atures (Tn) as well as the MI boundary above Tn. We 
present in §3.2 the results of calculations for the densities 
of states (DOS) as well as the MI transition below Tn- 
We will show that the dynamical effects tend to increase 
the critical Coulomb interaction for the MI transition. 
The numerical results for the ferromagnetic properties 
on the fee lattice in infinite dimensions are presented in 
§3.3. The Curie temperatures as well as the paramagnetic 
susceptibilities are shown to be described quantitatively 
by means of the dynamical CPA-I-HA in the quarter-filled 
regime of electron number. In the last §4 we summarize 
our numerical results. It is emphasized that the dynam- 
ical CPA is applicable to the quantitative description of 
the magnetic properties in high dimensions. 

2. Dynamical CPA in the ferro- and antiferro- 
magnetic states 

We adopt in the present paper the Hubbard model as 
follows. 

if = X] ^oriia + ^ UjAa'^jc + ^ ni^^i^ ■ (1) 

Here eo is the atomic energy level, tij is the transfer in- 
tegrals between sites i and j, and U is the intraatomic 
Coulomb interaction energy parameter. a\^{aiu) denotes 
the creation (annihilation) operator for an electron with 
spin (T on site i, and Uia = al^aja is the number operator 
for the electron. 

In the dynamical CPA,^^' we transform the two-body 



interaction Uni-^rii^^ in the free energy into a dynamical 
one-body potential Vi with time-dependent random spin 
and charge fields using the Hubbard-Stratonovich trans- 
formation. Then introducing the site-diagonal effec- 
tive potential Ej {i.e., the coherent potential) into the 
potential part of the free energy, we expand the correc- 
tion terms containing 7;^ — S; with respect to the site. The 
zeroth term is the coherent term which does not de- 
pend on the dynamical potential at all. The next term 
consists of the single-site terms each of which contains 
the dynamical potential on the same site. The higher- 
order terms A describe the inter-site spin and charge 
fluctuations. 

In the dynamical CPA, we adopt the single-site ap- 
proximation which neglects the higher-order inter-site 
corrections AT. The free energy of the dynamical CPA 
is then written as follows. 

^CPA = m - ^-'E'-J ^'^^^ ^-^""''^^'^ ■ (2) 

Here /3 denotes the inverse temperature, is the static 
exchange field on site i, and E^^^ (^j) is a single-site effec- 
tive potential projected onto the static field ^j. 

The effective potential E^^\^) consists of the static 
part E^l\^) and the dynamical part -E'dyn(0- 

^«(c) = £«(e) + <„(e). (3) 

The static potential E^l' (^) is obtained by neglecting the 
time-dependence of the field variables and is given by 

-\u{n,{0'-e), (4) 

vl^\^) = eo - M + luhii^) - hjia. (5) 

Here vf^{i^) is the Hartree-Fock type static potential, /z 
denotes the chemical potential. Yii„(i(jJi) is the frequency 
representation of a time-dependent coherent potential 
Sicr(T), T being the imaginary time, = (2/ -|- l)7r//3 
denotes the Matsubara frequency. The electron number 
fli(^) for a given exchange field ^ is defined by 

nm = -^Y.^iS^ui). (6) 

The Green function Gicr{iu)i,^) will be determined later 
self-consistently (see eq. (16)). 

The coherent Green function Fi„{icoi) in eq. (4) is de- 
fined by 

E,^{iuji) = [{iuJi -Ha~ S,(«w())"'] „: ■ (7) 

Here [Hq)^ = {cq — ^)5ij+tij{\ — 6ij) is the one-electron 
Hamiltonian matrix element in eq. (1) and = 

The dynamical potential E^^^{^) in eq. (3) is given in 
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the harmonic approximation (HA) as 
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Here gia{l) is the static Green function on site i, which 
is defined by 

~9ia{l) = [Fi^iiioi)-' - vl'^\0 + ^iaiicoi)] . (12) 



For the higher-order terms of B^^^^ , one can adopt the 
following form of the asymptotic approximation. 



E 



.fc=0 



(13) 



fe=0 



m=0 ^ 

^L°l(^,fc)CL(-^,fc) 

27r(^ — m)i/ _ 



+ 



z/3 



xb^^Uiy,k)b[]^^_,^{-i.,k)] . (14) 

The functions bmlri^i^, k) and 5mCT(±t^, k) in the 0-th or- 
der and the second order asymptotic approximation are 
given in Appendix B of Ref.^^^ 

The coherent potential 'Ei^{iuji) is determined so that 
the higher-order inter-site corrections become minimum. 
The condition called the CPA equation is given by 



(15) 



Here the average ( )eff at the l.h.s. (left-hand-side) of the 
above equation means taking a classical average with re- 
spect to the effective potential i?*^*'(^). The on-site dy- 



namical impurity Green function is given by 



EE 



Gi^{iujui) = gia{l) + - — is^s^ ■ (16) 

1+EE^i:^ 

Here Kia{iuji) = 1 - Fic,{iu>i)~'^SFic{iu]i)/6'Ei^{iuJi). 

The average electron number (rij) and the local mag- 
netic moment (mj) are given by 



{rii) = {ni{^i))ef[ 



(17) 



{■nii) = (ei)eff ■ (18) 

The double occupation number {rii^nn) and the ampli- 
tudes of local moment on site i are obtained from the 
following expressions. 

2 



dU 



(19) 



J f / eff 



(m^) = {m) - 2{ni^nii) . (20) 

Here [ ]„ means to take the derivative of a quantity fixing 
the static potential vl'^\^i). 

Equations (15) and (17) form a set of self-consistent 
equations to determine the coherent potential {'Eicr{iu)i)} 
and the chemical potential eo — /x. However it is time con- 
suming to solve the CPA equation (15) because one has 
to take the average ( )cn at each frequency uji. The fol- 
lowing decoupling approximation simplifies the numeri- 
cal calculations. 



9=±1 



Gicr{iwi,qxi) = Fia{iwi) . (21) 



Here Xi = (Cf )eff- Note that the approximation is cor- 
rect up to the second moment (i.e., (Ci)cff and (^j?)off). 

We assume in this work that each site is crystallo- 
graphically equivalent to each other. In the para- and 
fcrro-magnetic states, we can assume a sitc-indcpcndcnt 
coherent potential E(^(iwj). The site- independent coher- 
ent Green function F„(jbbJi) is then given by 

p(e)de 



Fa{l^l) = [ 



• (22) 

eo - ^a{iuJi) - e 

Here p(e) is the density of states (DOS) for the nonin- 
teracting Hamiltonian matrix tij. 

In the antiferromagnetic state with sublattice mag- 
netization, wc have two types of coherent potential; 
T,^\iuji) on the up-spin sublattice and \i'jJi) on the 
down-spin sublattice. Because of the symmetric relation 
SCT~^(iw;) = S^^(iw;), it is enough to require the self- 
consistency of T,'"c^\iLUi) on the up-spin sublattice. The 
corresponding coherent Green function Fi'^\iuii) is given 
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Fig. 1. Sublatticc iiiagnctizatioii vs. temperature curves and the 
amplitude of local moments {m^)^/'^ for the half-filled Hubbard 
model on the Bethe lattice in infinite dimensions. Coulomb inter- 
action energy parameter is chosen to be (7 = 0.5 in unit of half 
the band width for noninteracting density of states. Dashed curves 
show the results of the static approximation (SA), solid curves show 
the results of the harmonic approximation (HA), closed circle with 
error bar indicates the Neel temperature obtained by the Quantum 
Monte-Carlo (QMC) method.^^. 30) 
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(23) 



We adopted eqs. (22) and (23) in the para- and 
antifcrro-magnctic calculations on the Bcthc lattice and 
in the ferromagnetic calculations on the fee lattice in 
infinite dimensions. Moreover in the self-consistent cal- 
culations with use of the HA, we adopted eqs. (10) and 
(11) for and terms and took into account the 
higher-order terms up to U^^ using the asymptotic ap- 
proximation, i.e. eqs. (13) and (14). 

3. Numerical Results 

3.1 Antiferrom,agnetic states on the Bethe lattice 

We present in this subsection the numerical results of 
calculations for the antiferromagnetic states at half filling 
on the Bethe lattice in infinite dimensions. The DOS for 
the noninteracting system is given by a semi-elliptical 
form as 

= ^^/^^^^ • (24) 

Here 2W denotes the band width. In the present and 
next subsections, we adopt the energy unit to be W = 1. 

We found that the antiferromagnetism is stabilized at 
half filling with decreasing temperatures. The numeri- 
cal results of sublattice magnetization vs. temperature 
curves are presented in Fig. 1 for rather small Coulomb 
interaction energy parameter U = 0.5. Calculated Neel 
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Fig. 2. Sublattice magnetization vs. temperature curves and the 
amplitude of local moments for U = 3.5. Notations are the same 
as in Fig. 1. 



temperatures Tn are 0.017 in the static approximation 
(SA) and 0.011 in the harmonic approximation (HA), 
respectively. The dynamical effects reduce the ground- 
state sublattice magnetization (extrapolated value) by 
15% and Tn by 36%. The calculated Neel temperature 
with dynamical corrections seems to be in good agree- 
ment with the QMC result as shown in Fig. 1. Although 
the agreement might not be so convincing because of 
a large error bar in the QMC calculations, quantitative 
description in the weak Coulomb interaction regime is 
expected because the HA is exact up to the second order 
in U in the weak Coulomb interaction limit. The ampli- 
tude of local moment (rn^)^/^ is also enhanced by the 
dynamical effects because electron correlations suppress 
the double occupancy of electrons even at finite temper- 
atures. 

The results for a strong Coulomb interaction U = 3.5 
are presented in Fig. 2. The sublattice magnetization and 
Neel temperature are much enhanced as compared with 
the case oi U — 0.5. The ground-state sublattice magne- 
tization (0.960 12b for the SA and 0.955 hb in the HA) 
is close to the atomic value 1.0 /xb- Calculated Tn are 
0.0645 in the SA and 0.0654 in the HA, respectively; the 
Tn is slightly increased by the dynamical effects, but is 
smaller than the QMC value^^^ by 18%. 

We have calculated the Neel temperatures for various 
Coulomb interaction energy parameters. The results are 
summarized in Fig. 3. It should be noted that we adopted 
the decoupling approximation (21) in the calculations of 
Figs. 1 and 2 in order to reduce the computation time. To 
clarify the difference between the decoupling approxima- 
tion (21) and the full self-consistent CPA (15), we per- 
formed the full static calculations with use of (15). The 
calculated Tn vs. U curve is shown by solid curve and 
is compared with those in the decoupling approximation 
as well as the QMC results. In both the weak and strong 
Coulomb interaction regimes {i.e., U < 1 and U > 3), 
calculated Tn with different approximation schemes yield 
basically the same result. In the intermediate Coulomb 
interaction regime {i.e., 1 < U < 3), the Neel tempera- 
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Fig. 3. Nccl temperatures (Tn) in the SA with decoupling ap- 
proximation (dotted curve), the full SA (solid curve), the HA 
with dynamical corrections (closed squares), and the QMC method 
(closed circles).^"' The open squares around U = 3.0 above Tn 
indicate the metal-insulator crossover points in the full SA. Below 
Tn, the antiferromagnetic state (AF) is stabilized, while above Tn 
the paramagnetic metal (PM) and the paramagnetic insulator (PI) 
are realized. 



ttires in the decoupling approximation deviate from the 
full static results by several percent. It should be noted 
that Tn calculated by the full SA quantitatively agree 
with the QMC results. This indicates that the SA pro- 
vides us with a good starting point to calculate quantita- 
tively the magnetic properties of electrons especially at 
high temperatures. 

We have estimated the Neel temperatures in the full 
HA by adding ATn = Tn(HA) - TnCSA) in the decou- 
pling approximation to Tn in the full static one.^'') The 
results are shown in Fig. 3 by closed squares. We find 
that the dynamical effects in the harmonic approxima- 
tion reduce Tn in the weak Coulomb interaction regime 
{U < 1) and slightly increase Tn in the strong Coulomb 
interaction regime {U > 2). Calculated Tn are in agree- 
ment with the QMC results in the weak U regime, but 
are midcrcstimatcd in the intermediate and strong U 
regimes (1 < U). In particular the present dynamical 
calculations hardly correct Tn in the SA in the strong U 
regime (2 < U). One of the possible reasons for insuffi- 
cient corrections is that the decoupling approximation is 
not suitable at high temperatures; it tends to overesti- 
mate the magnetic entropy at high temperatures because 
the approximation implies to replace a broad distribution 
p(0 = exp(-^£;(')(0)/ / exp(-/3.B(*)(^)) with a two- 
delta function. The higher order dynamical corrections 
should also be taken into account in eq. (9) for strong U 
regime in order to make more reasonable corrections. 

3.2 Metal-insulator transition 

Single-particle excitation spectra are obtained by 
means of a numerical analytic continuation'^®^ of the self- 
energy for temperature Green function, i.e., the coherent 
potential T,i„{iuJi). It should be noted that the coherent 
potential obtained by solving the CPA equation (21) in 
the decoupling approximation is an approximate solution 



Fig. 4. Density of states (DOS) for U = 1.7 and T = 0.06 in 
the SA (dashed curve), and in the HA with dynamical corrections 
(solid curve). 
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Fig. 5. DOS for (7 = 2.2 and T = 0.06 in the SA (dashed curve), 
and in the HA with dynamical corrections (solid curve). 



to the full equation (15). In order to obtain accurately 
the DOS as the single-particle excitations, we improved 
the coherent potential adopting the average f-matrix ap- 
proximation (ATA)^'^^' after we solved eq. (21). 



E^^^(iw„) = Ei<,(iw„)-h 



{Gia {iC0„ ,^))eS-Fia- {iuj„ ) 



(25) 



{Gia- {iLOn,0)esFia{iiOn) 

Here the coherent potential in the decoupling approxima- 
tion was used at the r.h.s. (right-hand-side) of the above 
equation, but the full average ( )eflF of the impurity Green 
function is taken. The ATA is a one-shot correction to 
the fuU CPA (15). 

Making use of the Fade numerical analytic continua- 
tion,^®^ we obtained '£f'J'-^{uj + iS). Here uj is an energy 
variable on the real axis and i5 is an infinitesimal positive 
number. Using the self-energy, we obtained the coherent 
Green function Fa^\uj + i6), and calculated the DOS via 
the relation, 



p,{w) = --Im F^+\uj + i6) 



(26) 
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Fig. 6. Double occupation numbers as a function of U at T = 
0.02 in the Hartrcc-Fock (HF) approximation (dashed curve), the 
SA (dotted curve), and the HA (sohd curve). The paramagnetic 
state is assumed. Open squares show the double occupation num- 
bers in the QMC calculations^®^ at T = 0.03125. The vertical hnes 
indicate the critical Coulomb interaction Uc in the SA {Uc = 1.5: 
thin dotted line) and HA {Uc = 1.6: dot-dashed line). Both values 
of Uc are obtained within the decoupling approximation (21). 

The system changes from metal to insulator with the 
formation of a gap on the Fermi level with increasing 
Coulomb interaction. In the full SA, the gap is gradu- 
ally formed above the Neel temperatures. The crossover 
points are shown in Fig. 3 by open squares. When we as- 
sume the paramagnetic state below Tn, the formation of 
a gap becomes clearer at lower temperatures, indicating 
the metal-insulator (MI) transition. 

Figure 4 shows the DOS in the metallic state near 
the MI transition. We find the upper and lower Hubbard 
bands. A gap is almost opened on the Fermi level in the 
SA, while the quasiparticlc peak remain in the dynamical 
calculations. When the Coulomb interaction energy U 
is increased, a gap is opened as shown in Fig. 5. The 
dynamical effects shift the spectral weight to the lower 
energy region. 

The metal-insulator transition is accompanied by a lo- 
calization of electrons. We show in Fig. 6 an example of 
double occupation number {n^ni) vs U curve at low tem- 
peratures. The Hartree-Fock approximation yields the 
constant value 0.25 irrespective of U because of no on-site 
correlations. The double occupation number in the SA 
reduces to the Hartree-Fock value in the weak Coulomb 
interaction limit. It rapidly decreases near the MI transi- 
tion, and gives the alloy-analogy results of the Hubbard 
theory^) in the strong U limit. The dynamical effects 
based on the HA reduce further the double occupancy in 
the weak U regime and enhance it in the strong U regime 
as shown in Fig. 6. 

There are no QMC calculations on the double occu- 
pation number at T = 0.02. We plotted in Fig. 6 the 
QMC results^^) at T = 0.031 for qualitative or semi- 
quantitative comparison with the present results. The 
QMC value at U — 0.25 indicates that the dynami- 
cal CPA-I-HA yields the quantitative result in the weak 



Fig. 7. Phase boundaries of metal-insulator transition below Tn 
in the SA (open squares), the HA (closed squares), and the QMC 
method (closed circles).'^-') The critical Coulomb interaction in the 
alloy-analogy approximation (AA)®^ at T = is indicated by the 
open circle, while the exact results of Ud and Uc2 at T = based 
on the projective self-consistent technique (PSCT)^^) are indicated 
by the stars. 

Coulomb interaction. The QMC result at U = 1.0 sug- 
gests that the HA is not sufficient to reduce the double 
occupancy in the intermediate regime of the Coulomb 
interaction. It should be noted that the critical Coulomb 
interaction in the QMC is Uc = 2.4 so that the data in 
2.0 < U < 2.4 are for metallic region. Thus, they are not 
directly compared with the present results in the insula- 
tor regime. The dynamical CPA-I-HA result at U = 2A 
agrees with the QMC result. These results suggest that 
the HA gives the quantitative results of (n^ni) in both 
the weak {U < 1.0) and strong {U > 2.4) Coulomb inter- 
action regimes. In the intermediate Coulomb interaction 
regime (1.0 < U < 2.3), we need to take into account 
more dynamical effects to obtain (n-i-nj,) quantitatively 
at low temperatures. 

We have plotted in Fig. 7 the MI transition points 
below Tn being obtained from the gap formation in 
the DOS. The critical Coulomb interactions Uc for gap 
formation in the full SA increase with increasing tem- 
perature, and merge into the QMC results^^^ around 
T « 0.06. Since the full SA is exact in the high temper- 
ature limit, this indicates that the full SA quantitatively 
describes Uc at T > 0.06. It should be noted that these 
data are smoothly connected to the data points above 
Tn for a gap formation in Fig. 3. Below T = 0.06, the Uc 
in the SA deviate from the QMC results, and reduces to 
Uc = 1.0 in the Hubbard alloy-analogy approximation^^ 
at T = 0. 

The dynamical results of Uc in Fig. 7 are obtained 

by adding the dynamical correction AUc = t/c(HA) — 
(7c (S A) in the decoupling approximation to C/c(SA) in 
the full SA.^^) It is known that there are two critical 
Coulomb interactions at T = 0, Ud for the gap forma- 
tion and Uc2 for the disappearance of the quasiparticlc 
state. But we did not find numerically two solutions 
near Uc at finite temperatures (T > 0.01). The critical 



J. Phys. Soc. Jpn. 



FULL PAPERS 



- 








- 






- 


- 












M 








□\ 






n=0.58 








U=4 










r-i 1/ 1 





0.02 0.04 0.06 0.08 0.1 0.12 



T 

Fig. 8. Magnetization (Af), inverse susceptibility (x^'') and am- 
plitude of loeal moment {l^vii^)^!^') as a function of temperature 
(T) on the fee lattiee in infinite dimensions. Solid curves: dynami- 
cal CPA+HA, dotted curves : dynamical CPA+SA. Open squares 
(open circles) are the magnetizations (inverse susceptibilities) in 
the 

Coulomb interactions in the HA are larger than those 
of the full static approximation. The difference between 
the two becomes smaller with increasing temperature 
and vanishes at T > 0.06 as it should be. The critical 
Coulomb interaction extrapolated to the zero temper- 
ature is about 1.6 in the present HA, while the exact 
values based on the projective self-consistent technique 
(PSCT)32) g^j.g reported to be C/d = 2.40 and C/c2 = 2.84. 
One has to take into account the dynamical terms higher 
than more seriously in order to describe the MI tran- 
sition quantitatively. 

3.3 Ferromagnetism on the fee lattiee in infinite dimen- 
sions 

We have investigated the ferromagnetic properties of 
the fee lattice in infinite dimensions by means of the dy- 
namical CPA-I-HA in order to examine the quantitative 
aspects of the theory. The fee lattice in d dimensions is 
defined as a lattice with 2d{d — 1) nearest neighbors ex- 
pressed hy R = ±ei± Sj with two different cubic unit 
vectors and Bj = 1, 2, • • • ,d). In infinite dimen- 
sions the noninteracting DOS is given by^^^ 

(1 + V2e) 

p{e) = \ ^ . (27) 
^7r(l + V2e) 

Here the energy unit is chosen so that the bottom of band 
edge is given by — 1/\/2. 

Note that the DOS (27) monotonically increases with 
decreasing energy e and diverges at the band edge e = 
— l/y/2. Therefore the DOS is expected to be favorable 
for the ferromagnetism in the low density region. Since 
the numerator in the DOS is reduced by a factor of 1/e 
at e = 1/ -s/2, we may regard W = \f2 as a characteristic 
band width of this system. The coherent Green function 
Fcr{iu!i) was calculated from eq. (22) via numerical inte- 
gration with respect to energy e. 



Figure 8 shows an example of calculated magnetiza- 
tion vs temperature curves in various approximations 
for electron number n = 0.58 and Coulomb interaction 
U = 4.0. Extrapolated ground-state magnetizations are 
almost saturated in both calculations. The magnetiza- 
tion vs. temperature curves considerably deviate down- 
wards from the Brillouin curve. The SA overestimates 
the Curie temperature (Tq = 0.0800) as compared with 
the QMC^^) (Tc = 0.0512). The dynamical effects based 
on the HA reduces Tq by 30% and rcsuhs in Tc = 0.0559. 
The result agrees with the QMC^^) (Tc = 0.0512) within 
9% error. We have also calculated the paramagnetic sus- 
ceptibilities above Tc by adding the infinitesimal mag- 
netic field. Calculated susceptibilities follow the Curie- 
Weiss law; X = '^eff/(-^ ~ ^c)- Here rries is a constant 
called effective Bohr magneton number. Calculated effec- 
tive Bohr magneton numbers are 1.00 /zb in the SA, 1.13 
fXB in the HA, and 1.18 /xb in the QMC, respectively. We 
find a good agreement between the dynamical CPA-I-HA 
results and the QMC ones. 

We have calculated Tc as a function of electron num- 
ber n by extrapolating the inverse susceptibilities to 
lower temperatures at various n. The results are pre- 
sented in Figs. 9 and 10 for an intermediate Coulomb 
interaction strength U = 2 and a strong Coulomb in- 
teraction U = 4, respectively. In the SA, we obtain a 
finite value of Tc for infinitesimal electron number n. 
The Tc increases linearly with increasing n, and shows 
a maximum Tc = 0.053 (0.083) at n^ax = 0.60 (0.70) 
fov U = 2 {U = 4). When electron number n ap- 
proaches to one, Tc vanishes again at n = 0.93 (0.98) 
for t/ = 2 {U = 4). The dynamical effects reduce Tc, 
for example, by 35 % (25 %) at n^ax = 0.60 (0.70) 
for U = 2 (U = A). As shown in Figs. 9 and 10, the 
dynamical results quantitatively agree with the QMC 
results around the quarter filled electron number {i.e., 
0-4 < n < 0.6) for both U = 2 and t/ = 4. However, 
the Tc vs n curves shift to the higher density region by 
about An = 0.1; calculated Tc in the low density regime 
(0 < n < 0.3)) arc tmderestimated, and those in the 
half-filled regime (0.7 < n < 1) are overestimated. 

The behavior of Tc for small n is sensitive to the 
approximation. In the Hartree-Fock approximation, Tc 
should be finite for infinitesimal number of electron ac- 
cording to the Stoner condition p{0)U > 1 because the 
noninteracting DOS diverge at n = 0. This holds true 
even in the static approximation because the SA reduces 
to the Hartree-Fock approximation in the weak Coulomb 
interaction limit. The QMC results suggest that Tc van- 
ish at finite electron number in the low density region. 
There is no exact analysis on the low-density behavior of 
Tc for this system as far as we know. The present calcu- 
lations indicate that Tc vanishes at n = 0.2 for [/ = 2, 
and at n = 0.3 for {7 = 4. 

Magnetic properties at the high temperatures above 
Tc arc expected to be described well by the dynamical 
CPA-I-HA. We present in Figs. 11 and 12 calculated cf- 
fectivc^ Bohr magneton number (mofr) and amplitude of 
local moment ((m^)"'^/^) as a function of electron number 
n. Note that {m^y/^ = V3{m^y/'^ where (m2)V2 jg the 
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Fig. 9. Curie temperature (Tc) vs. electron number (n) at [/ = 2 
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Fig. 11. Calculated effective Bohr magneton numbers (meff) as 
a function of the electron number n in the HA (closed circles) 
and in the SA (open circles). Coulomb interaction energy is fixed 
to be [/ = 2. The amplitudes of local moments {m^)^^'^ axe also 
presented by open squares. 
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Fig. 10. Curie temperature vs. electron number at 1/ = 4 in the ri 
SA (dotted curve), the HA (closed circles), and the QMC (open 

squares with error bars ).33' Fig. 12. Same as in Fig. 11 but for U = 4. The QMC result at 

n = 0.58 is shown by +. 



amplitude of local moment for z component given by eq. 
(20). These quantities were obtained at high tempera- 
tures (T w 0.09). In the case ofU = 2, the effective Bohr 
magneton number monotonically increases with increas- 
ing electron number from to 1, and shows the maxi- 
mum 1.64 hb at half-filling. Although rricff agrees with 
the amplitude (m^)^/^ in the atomic limit, it is in gen- 
eral smaller than the latter in the metallic state. At half- 
filling, the difference is small; nicg « (m^)^/^ because of 
considerable suppression of charge fluctuations. By com- 
paring the effective Bohr magneton numbers in the HA 
with those in the SA, we find that the dynamical effects 
enhance the effective Bohr magneton number, though the 
SA describes nies quantitatively at half-filling. 

For larger value of U = 4, we find similar behavior of 
TOeff as a function of n, but the effective Bohr magneton 
number shows a maximum at n « 0.95, and jumps from 
1.3 fiB to 1.6 /^B at n = 1. Sudden jump to the value 1.6 
/xb being close to the atomic value a/S is associated with 
the formation of an insulator at half-filling. Calculated 
effective Bohr magneton number at n = 0.58 quantita- 



tively agrees with the QMC result as shown in Fig. 12. 

4. SumrriEiry 

In the present paper we have investigated magnetic 

and electronic properties of the Hubbard model in infi- 
nite dimensions on the basis of the dynamical CPA com- 
bined with the harmonic approximation (HA) in order to 
clarify the accuracy of the theory at finite temperatures. 

In the dynamical CPA, we transform the interacting 
electron system into an independent electron system cou- 
pled with the timc-dcpcndcnt random charge and ex- 
change fictitious fields on the basis of the functional in- 
tegral method. Introducing the coherent potential, we 
made a single-site approximation in the free energy. The 
coherent potential is self-consistently determined by the 
condition that the average of dynamical impurity Green 
function embedded in the effective medium (i.e., the co- 
herent potential) should be equal to the coherent Green 
function for the uniform medium. The former Green 
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function is calculated by means of the HA. The dynam- 
ical CPA+HA becomes exact in the high-temperature 
limit, leads to the free energy being exact up to the sec- 
ond order in the weak U limit, and yields the exact result 
in the strong U limit in infinite dimensions. 

We have performed the numerical calculations with 
use of the dynamical CPA+HA in which the dynami- 
cal corrections have been taken into account up to the 
4th order in U exactly and the higher-order corrections 
up to the 16th order in U have been taken into account 
with use of the second order asymptotic approximation. 
In the present calculations for the half-filled band Hub- 
bard model on the Bcthc lattice in infinite dimensions we 
found that the Neel temperatures Tn are well described 
by the zeroth approximation to the dynamical CPA, i.e., 
the static approximation (SA). The dynamical correc- 
tions based on the HA reduce further Tn in the weak 
Coulomb interaction regime, leading to a quantitative 
agreement with the QMC results. In the intermediate 
Coulomb interaction regime, the dynamical corrections 
in the present approximation tend to underestimate Tn 
by about 10 %. In the strong Coulomb interaction regime 
we found that the dynamical corrections slightly increase 
Tn- These results suggest that the dynamical CPA+HA 
can describe Tn within 10% error in infinite dimensions. 

We have also verified using the same model at half 
filling that the SA quantitatively describes the MI 
crossover above Tn- It however underestimates the crit- 
ical Coulomb interaction Uc at low temperatures. The 
dynamical corrections based on the HA enhance Uc. But 
the corrections in the present version are not enough to 
reproduce quantitatively the critical Uc obtained by the 
QMC; the present theory leads to Uc smaller than those 
of the QMC by 30% at low temperatures (T - 0.02). 

We have investigated the ferromagnetic properties on 
the fee lattice in infinite dimensions. In this case, the 
SA overestimates the Curie temperatures Tq more than 
30% irrespective of electron number. The dynamical 
CPA+HA reduces Tc in the SA, and quantitatively de- 
scribes the Curie temperatures as well as the susceptibili- 
ties in the quarter-filled regime (0.4 ^ n < 0.6). However 
the Tc vs n curves tend to shift to the higher density re- 
gion by about An = 0.1 as compared with those in the 
QMC calculations, so that Tq in the low density region 
(0 < n < 0.3) are somewhat underestimated and those 
in the half-filled regime (0.7 ^ n < 1) are overestimated. 
We found that the effective Bohr magneton number rrics 
is smaller than the amplitude (m^)^/^ in itinerant elec- 
tron system except at half-filling where electrons tend to 
localize due to electron correlations for the intermediate 
Coulomb interaction strength. The dynamical effects en- 
hance nics for non half-filled bands, while the corrections 
are negligible at half-filling. 

We have to perform the full CPA calculations with- 
out decoupling approximation (21) in order to obtain 
more solid conclusions. Nevertheless, the present results 
of calculations indicate that the dynamical CPA com- 
bined with the HA is suitable for the quantitative cal- 
culations of magnetic properties at finite temperatures 
as well as the MI transition at high temperature region. 



It is not enough for the quantitative description of the 
MI transition at low temperatures and for the quantita- 
tive description of the ferromagnetism in the low-density 
regime as well as in the half-filled regime. For more quan- 
titative description in the intermediate Coulomb interac- 
tion and low-temperature regime, one has to treat more 
seriously the higher-order dynamical corrections. The im- 
provement is left for future work. 
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